g-breathers in finite two- and three-dimensional nonlinear acoustic lattices 
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Nonlinear interaction between normal modes dramatically affects energy equipartition, heat con- 
duction and other fundamental processes in extended systems. In their celebrated experiment Fermi, 
Pasta and Ulam (FPU, 1955) observed that in simple one-dimensional nonlinear atomic chains the 
energy must not always be equally shared among the modes. Recently, it was shown that exact 
and stable time-periodic orbits, coined g- breathers (QBs), localize the mode energy in normal mode 
space in an exponential way, and account for many aspects of the FPU problem. Here we take the 
problem into more physically important cases of two- and three-dimensional acoustic lattices to find 
existence and principally different features of QBs. By use of perturbation theory and numerical 
calculations we obtain that the localization and stability of QBs is enhanced with increasing system 
size in higher lattice dimensions opposite to their one-dimensional analogues. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 



Nonlinearity induced interaction between normal 
modes of extended systems is crucial for many funda- 
mental dynamical and statistical phenomena like ther- 
malization, thermal expansion of solids or turbulence in 
liquids. It is also important in many artificial systems 
where one aims at controlling the energy flow among 
the normal modes, preventing or efficiently channeling 
energy pumping due to resonances with external per- 
turbations etc. Among the many accumulated results 
in this area a seminal one is due to Fermi, Pasta and 
Ulam (FPU) who reported on the absence of thermaliza- 
tion of chains of atoms connected by weakly nonlinear 
springs jj|. They observed that the energy of an ini- 
tially excited normal mode with frequency ujq and wave 
number q did not spread over all other normal modes, 
staying almost completely locked within a few neigh- 
bouring modes in the normal mode space [3, y| . Longer 
waiting times yielded another puzzle of energy reccur- 
rence to the originally excited mode. Many efforts to 
explain the FPU paradox resulted in an extraordinary 
progress in this field: the observation of solitons [j|, 
size-dependent stochasticity thresholds pj, nonlinear res- 
onancesldL KAM tori Arnold diffusion, and many other 
issues Uli, 0, E3, El- The efforts to carry all these 
concepts into two-dimensional lattices have been also re- 
ported [13. However, according to Ford |2|, despite the 
richness of new topics, the original FPU problem was still 
waiting to be understood. 

Recently it was shown that the major ingredients of the 
FPU problem can be addressed within the promising con- 
cept of g-breathers (QBs), which are exact time-periodic 
solutions in the nonlinear FPU chain |l3|. These solu- 
tions are exponentially localized in the q-space of normal 
modes and preserve stability for small enough nonlin- 
earity. They continue from their trivial counterparts for 
zero nonlinearity at finite energy. In that limit they cor- 



respond to one mode with the seed wave number go being 
excited, and all the other modes being at rest. The sta- 
bility threshold of QB solutions coincides with the weak 
chaos threshold in [3. Persistence of exact stable QB 
modes surrounded by almost quasiperiodic trajectories 
explains the origin of FPU recurrences and the absence 
of energy equipartition. The scaling of the localization 
exponents with the relevant control parameters relates 
to results on higher order nonlinear resonance overlaps 
[3|. But perhaps the most important result was, that it 
needs only one ingredient to obtain QBs in FPU chains: 
a discrete nonequidistant frequency spectrum of normal 
modes, as induced by a finite system. Thus the door is 
opened to apply the concept of QBs to a variety of non- 
linear finite systems - truly common objects in nature 
and applications. One of such challenges is an extension 
of the notion of QBs into more physically realistic two- 
and three-dimensional FPU-type systems. 

In this Letter we report on the existence and remark- 
able properties of g-breathers in finite two-and three- 
dimensional nonlinear acoustic lattices. For fixed energy, 
nonlinearity coefficient and seed wave number the QB 
localization length stays finite in the 2d case and tends 
to zero in the 3d case with increasing size of the system. 
For both 2d and 3d cases the nonlinearity coefficient at 
the QB stability threshold increases in the same limit. 
This comprises a crucial difference from the Id case |l.l | 
where QBs delocalize and the nonlinearity threshold of 
the instability tends to zero in the limit of large chains. 

We consider quadratic and cubic lattices of N'^ (d — 2 
and 3 respectively) equal masses coupled by nonlinear 
springs with the Hamiltonian 



H = lT.ip'^ 
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where Xn{t) is the displacement of the n = (ni, 
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th particle from its original position, Pnit) its momen- 
tum, D(n) is the set of its nearest neighbors, and fixed 
(zero) boundary conditions are taken: Xn = ii ni — 
OT ni ^ N + 1 for any of the components of n. 
A canonical transformation 
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takes into the reciprocal wave number space with N^ 
normal mode coordinates Qq{t) = Qqi,...,qa{t)- The nor- 
mal mode space is spanned by q and represents a d- 
dimensional lattice similar to the situation in real space. 
The equations of motion then read 



^iQq = - 
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Here ^q = 4 ^ cj^. are the squared normal mode fre- 

1=1 
quencies with cj^. = sin (7rqi/2(A^ -I- 1)). Note, that all 

linear modes but the diagonal ones Qq^=,„=q^ are at least 
d-fold degenerate with respect to their frequencies. The 
coupling coefficients Cq^p^r,s induce a selective interac- 
tion between distant modes in the normal mode space 
similar to the Id case. 

For small amplitude excitations the nonlinear terms in 
the equations of motion can be neglected, and according 
to (EJ the g-oscillators get decoupled, each conserving its 

harmonic energy Eg = ^ ( Q^ + ^gQg j in time. Espe- 
cially, we may consider the excitation of only one of the 
g-oscillators, i.e. Eg j^ ior q = q^ only. Such exci- 
tations are trivial time-periodic and q-localized solutions 
(QBs) for /3 = 0. 

For the Id chain, such periodic orbits can be contin- 
ued into the nonlinear case at fixed total energy |lj| be- 
cause the nonresonance condition nQ.q^ ^ ^q#qo (here 
n is an integer) holds for any finite size |lj|. This ar- 
gument can be used straightforwardly for d = 2, 3 and 
large seed wave numbers q on the main diagonal such 
that 2Slq > 4. For all other seed wave numbers on the 
main diagonal we checked that the nonresonance condi- 
tion holds as well. For seed wave numbers off the main 
diagonal the above mentioned d-fold degeneracies could 
be lifted by considering anisotropic lattices. In fact these 
degeneracies are lifted by the nonintegrability of the non- 
linear lattice (EJ , and a discrete set of periodic orbits will 
remain. QBs are successfully observed in numerical ex- 
periments in the presence of the degeneracy and we did 
not find substantial difficulties in computing them. 

In the following we compute QBs as well as their Flo- 
quet spectrum numericall y us ing the same algorithms as 
for the Id atomic chain [lj|, and compare the results 
with analytical predictions, derived by means of asymp- 
totic expansions. 




FIG. 1: A QB for d = 2 with q^ = (3, 3) for (3 = 0.5, Etot 
1.5, frequency tl ^ 0.403, N = 32. 




FIG. 2: The normal mode energies along (a) the diagonal 
q{n) = nq^ and (b) the side direction p{n) = {qo,i,nqo,2) for 
a QB with qp = (3, 3), /3 = 0.5, Etot = 1.5, N = 16, 32, 48 and 
the analytical estimate Q (dashed lines). 



First, we turn to the case d = 2. We obtain various 
symmetric (with (go)i = (90)2, Fig lll2ll and asymmetric 
(with {qo)i 7^ (90)2) QBs in the lower frequency range, 
which are exponentially locaHzed in q-space. Note, that 
due to the parity symmetry of the model (Eq. Q is invari- 
ant under Xn -^ —Xn for all n) only the modes with odd 
components (91,(72) are excited by the Qq = (3i3) mode. 
In contrast to d = 1, the decay of the energy distribution 
(especially along the diagonal q{n) = (2n— l)qo) remains 
almost constant with increase of the lattice size (FiglSJ . 

The abovementioned degeneracy of the frequency spec- 
trum supports the existence of multi-mode QBs, namely 
those, which have two (or more) excited seed modes. 
Indeed, we continued multi-mode periodic solutions of 
the linear lattice Eg ^ for q G 5'(qo) = {? ■ ^g = 
riq^} into the nonlinear regime. For example, the set 
Qg = (2, 3), Qq = (3, 2) allows for two (asymmetric single- 
mode) QB solutions with the energy mainly concentrated 



in one of the two seed modes and two symmetric multi- 
mode QB solutions with the same energy in each of the 
seed modes, and oscillations being in- or out-of-phase. 

By an asymptotic expansion of the solution to (pj in 
powers of the small parameter a = I3/{N+1)'^ (in analogy 
to [13I) we estimate the shape of a QB solution Qq{t) with 
a low-frequency seed mode number Qq. The energies of 
the modes on the diagonal of the QB Qq, ^Qq,. . . ,(2n -|- 
l)qov< iN,N) read 
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Dashed lines in FigEta) are obtained using I^J and show 
very good agreement with the numerical results. The en- 
ergy distribution between other modes, involved in the 
QB is more compHcated, but the decay along the diag- 
onal Q gives a good estimate for it (Figl^b)). Note, 
that the localization along the diagonal is the weakest, 
at least for large N . The shape of the QB in the g-space 
then possesses the following properties: (i) the localiza- 
tion remains constant when the lattice size tends to infin- 
ity with all other parameters fixed, in contrast to the Id 
case where QBs delocalize as A oc (A^-f 1), (ii) in the limit 
of constant energy density e = Eq^/{N + 1)^ and wave 
vector of the QB Kq = Qo/i^ + 1) the localization re- 
mains constant for A^ ^ 00 (similar to the Id case), (iii) 
for smaller /3, Eq^ and larger Qq QBs compactify. These 
results are a consequence of a more fundamental scaling 
property of finite to infinite systems. 

We analyze the stability of the obtained periodic or- 
bits with standard methods of Hnearizing the phase space 
flow around the solutions and computing the eigenvalues 
and eigenvectors of the corresponding symplectic Floquet 
matrix [l3|. The orbits are stable if all complex eigen- 
values lie on the unit circle. The absolute values of the 
Floquet eigenvalues of QBs for symmetric Qq = (3, 3) and 
asymmetric Qq — (2, 3) are plotted versus (3 for different 
system sizes TV in FigEI Similar to the Id case, QBs 
are stable for sufficiently weak nonlinearities. When /3 
exceeds a certain threshold /?* , some eigenvalues get ab- 
solute values larger than unity (and some of them less 
than unity) and the QB becomes unstable. In remark- 
able contrast to the Id case, /?* rapidly increases with the 
size of the system FigEJa). For a series of computation- 
ally accessible large lattices N = 20,30,40 (not plotted 
in FiglSIa)) we found the QB with Qq = (3,3) to be 
stable at least up to /? = 10.0. For insufficiently large 
N the dependence f3*{N) may become non-monotonous 
(FigEtb)). It may be quite sensitive to the chosen seed 
wave number Qq, compare FigsE^a) and (b). 

The observed instabilities can be traced analytically, 
similar to the Id case. Using standard secular per- 
turbation techniques we approximate the frequency of 
the QB solution as Cl = f7q^(l -I- 9hp) + 0{h?), where 
h = 3f3Eq^/{N -t- 1)^ is a small parameter and p = 




0.995 



FIG. 3: Absolute values of Floquet eigenvalues (symbols) and 
analytical curves (solid lines) for QBs with Etot = 1.5 and (a) 
Qq = (3, 3),Af = 13,14,16 (the instabilities correspond to 
k = (1, 2) and k = (2, 1)) and (b) Qq = (2, 3), iV = 20, 25, 30 
(the instabilities correspond to fe = (1, 1) for TV = 20, 25 and 
fc = (1,2) for TV = 30). 
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Linearizing the equations of mo- 
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tion Q around a QB solution Qq = Qqit) + ^q we obtain 



i<l + ^l^<l = -^h{l + COS2nt)J2Cq,q„q„p/^l^p + 0{h^) 

P 

(5) 
The strongest instability, caused by primary parametric 
resonance in Q , comes from pairs of resonant modes q + 
p = 2qQ with a nonzero vector k — q — Qq — qo — p. The 
bifurcation point and the absolute values of the Floquet 
multipfiers involved in the resonance in its vicinity are 
represented by a complex expression, which demonstrates 
good agreement with the numerical results (FigEJ. We 
assume |A:i^2| << (90)1,2 << TV and approximate 



3i3*Eqj{N + 1)2 = sirtq + np- 2nqj/n. 
(90)2^1 - (90)1^2 
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This estimate explains the following basic features of the 
observed instability. The instability of the type fc that 
minimizes 1(90)2^1 — {qo)ik2\ is the first to occur. It re- 
sults in a non-monotonous and discontinuous dependence 
of /3*(qo, TV) for small lattices. It monotonously increases 
with TV while k is constant (FigE^a)); when k changes 
/?* changes as well (FigEJb)). Besides, the instability 
threshold scales as /3* ex q^^ and as /?* oc TV^. 

The 2d lattice supports also the existence of QBs with 
Qo located in the intermediate and high-frequency parts 
of the normal mode spectrum. 

Let us turn to the case d ~ 3. We again compute QB 
solutions as in the 2d case. A similar analysis shows, 
that for the low-frequency seed mode Qq the decay of 
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FIG. 4: The structure of the QB with Qq = (2, 2, 2) on the 
three-dimensional lattice A'^ = 11, /3 = 0.5, Btot = 1.5, fJ « 
0.897. The size of spheres is a linear function of the decimal 
logarithm of the linear mode energy Eg . 



the normal mode energies along the leading direction 
3qo,. . . ,(2n + l)^^,. . . < {N, N, N) is given by Q with 
d = 3, and fits well the shape of the numerically com- 
puted QBs with Qo = (3,3,3) (FigEJ. In contrast to 
the 2d case, the localization length even decreases with 
increasing lattice size X^d oc (A^ + 1)~^. For constant 
energy density £ = E/{N + 1)'^ and wave vector of the 
QB Kq, the degree of localization is independent of the 
lattice size, as it is for the Id and 2d cases. Despite hav- 
ing difficulties when calculating the stability of QBs in 
sufficiently large 3d lattices due to limited machine per- 
formance, the analytical treatment has been done similar 
to lower lattice dimensions. We find that (5* on {N + 1)^ 
for constant energy of the lattice and does not increase 
for fixed e and q^. We also predict its sensitive depen- 
dence upon the choice of the seed mode number of the 
QB and the size of the lattice when being small. 

The main reason for the crucial improvement of local- 
ization and stability of QBs in higher dimensions relies 
on the fact that the effective strength of nonlinear inter- 
mode coupling decreases exponentially due to the factor 
(iV+1)-'^. 

In conclusion, we report on the existence and remark- 
able properties of g-breathers as exact time-periodic so- 
lutions in nonlinear two- and three-dimensional acoustic 
systems, thus extending the concept of one-dimensional 
QBs to more physically relevant objects. They are expo- 
nentially localized in the g-space of the normal modes and 
preserve stability for small enough nonlinearity. In the 
limit of infinite lattice size the localization length stays 
constant for d = 2 and tends to zero for d = 3 when 
fixing Qo cind the total energy. At the same time the 
localization length does not depend on the lattice size 
N when fixing the energy density and the wave vector 
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FIG. 5: The normal mode energies along the leading direction 
q{n) = (n, n, n) for q^ — (3, 3, 3),Etot = 1.5, A*' = 20 and the 
analytical estimate (^J (dashed lines). 



Kq = Qo/iN + 1). In contrast to the one-dimensional 
case the observed instability threshold in the nonlinear 
coupling strength increases with increasing lattice size, 
and stays constant if the energy density is fixed. 

M.I., O.K., and K.M. appreciate the warm hospitality 
of the Max Planck Institute for the Physics of Complex 
Systems. M.I. and O.K. acknowledges support of the 
"Dynasty" foundation. 



[1] E. Fermi, J. Pasta, and S. Ulam, Los Alamos Report 

LA-1940, (1955); also in: Collected Papers of Enrico 

Fermi, ed. E. Segre, Vol. II (University of Chicago Press, 

1965) p. 977-978; Many-Body Problems, ed. D. C. Mattis 

(World Scientific, Singapore, 1993). 
[2] J. Ford, Phys. Rep. 213, 271 (1992). 
[3] Chaos 15 Nr.l (2005), Focus Issue The Fermi-Pasta- 

Ulam problem - The first fifty years, Eds. D. K. Campbell, 

P. Rosenau and G. M. Zaslavsky. 
[4] N. J. Zabusky and M. D. Kruskal, Phys. Rev. Lett. 15, 

240-243 (1965). 
[5] F. M. Izrailev and B. V. Chirikov, Institute of Nuclear 

Physics, Novosibirsk, USSR, 1965 (in Russian); Dokl. 

Akad. Nauk SSSR 166, 57 (1966) [Soviet. Phys. Dokl. 

11, 30 (1966)]. 
[6] B. V. Chirikov, Atomnaya Energia 6, 630 (1959) [Engl. 

Transl. J.Nucl. Energy Part C: Plasma Phys. 1, 253 

(I960)]. 
[7] J. De Luca, A. J. Lichtenberg, and M. A. Lieberman, 

CHAOS, 5, 283 (1995). 
[8] D.L. Shepelyansky, Nonlinearity 10, 1331 (1997). 
[9] P. Bocchierri, A. Scotti, B. Bearzi, and A. Loigner, Phys. 

Rev. A 2, 2013 (1970); L. Galgani and A. Scotti, Phys. 

Rev. Lett. 28, 1173 (1972); A. Patrascioiu, Phys. Rev. 

Lett. 50, 1879 (1983). 
[10] H. Kantz, Physica D 39, 322 (1989); H. Kantz, R. Livi 

and S. Ruffo, J. Stat. Phys. 76, 627 (1994). 
[11] L. Casetti, M. Cerruti-Sola, M. Pettini and E. G. D. Co- 
hen, Phys. Rev. E 55, 6566 (1997). 
[12] G. Benettin, Chaos 15, 015108 (2005) and references 

therein. 
[13] S. Flach, M. V. Ivanchenko and O. I. Kanakov, Phys. 

Rev. Lett. 95, 064102 (2005); S. Flach, M. V. Ivanchenko 

and O. I. Kanakov, Phys. Rev. E 73, 036618 (2006). 
[14] J. H. Conway and A. J. Jones, Acta Arith. XXX, 229 



(1976). 



